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We discuss testing improved actions in the context of finite volume gauge theories, 
where both results for the continuum and the Wilson lattice action are known 
analytically for volumes up to 0.7 fermi across. A new improved action is in- 
troduced, obtained by adding a 2 X 2 plaquette to the Liischer-Weisz Symanzik 
action, for which the gauge field propagator greatly simplifies. We call this the 
square Symanzik action. We present the tree-level parameters of this improved 
action and the value of its Lambda parameter. We also give some Monte Carlo 
results and discuss some of the issues related to violations of unitarity at the scale 
of the lattice cutoff due to next-to-nearest coupling in the time direction. 



1 Introduction 

Lattice gauge theory has proven to be a viable tool for non-pcrturbative 
studies in field theories and non-Ahelian gauge theories in particular. In per- 
turbation theory it can be provenlau to have a continuum limit, and all results 
indicate that non-perturbatively the same will holdoB. Despite the tremen- 
dous increase of computer power since the first Monte Carlo calculations 
it remains a technical challenge to make the lattice spacing a small and the 
physical volume large enough. In particular problems of critical slowing down 
at small values of the coupling, implying that the algorithm is not efficient 
in probing the relevant portions of field space, makes a straightforward ap- 
proach of decreasing the lattice spacing (while keeping the physical volume 
large enough) a costly procedure. Instead one can use the increased computer 
power of today to study alternative lattice actions, which are chosen to remove 
as much as possible the scaling violations introduced by a finite lattice cut- 
off. The computational overhead in using a more complicated lattice action is 
usually not more than a few to ten times the cost for the Wilson action. 

1.1 Symanzik improvement 

The notion of improved actions was introduced more than a decade ago and i: 
particular the study of Symanzik u for scalar theories and Liischer and Weisz 

"Talk presented by the last author at the second workshop "Continuous advances in QCD" 
Univ. of Minnesota, March 28-31, 1996. 
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for non-Abelian gauge theories have been influential. The last few years have 
seen a surge in applying these ideas in actual simulations!! When considering a 
lattice action, here restricting ourselves to pure gauge theories, the dynamical 
variables live on the (directed) links of the lattice and are elements of the gauge 
group Ufj, (x) . The connection to a continuum configuration is provided in terms 
of parallel transport of the vector potential (x) along the link x — > x + /t (/} 
the directional vectors of the hypercubic lattice) 

U^x) = Pexp( [ A^(x + sp,)ds) . (1) 
Jo 

This allows us to express any of the lattice actions in powers of the lattice 
spacing. Terms that are of order a 4 (in four dimensions), being the volume 
of a unit cell on the lattice, correspond to the naive continuum limit. Higher 
order terms of which we will quote the powers of a relative to this volume 
factor, correspond ta.irrelevant operators. For example, for the Wilson action 
one finds the result lid 

S = ^Tr(l-f/(P))= ^^(l-Q) 

P x,n,u xt'Jl 

= Tr (! - U»{x)Uu{x + fi)Ul(x + v)Ui{x)) 

X,(J,,l> 

+ ^ {(PlF^x)) 2 + (VlF^(x)) 2 ) + X -VlF^{x)VlF^{x) 

-\{V,V u F,„{x)f}]+0{a w ) , (2) 

Note that apart from the corrections to the action density, sums need to be 
converted to integrals (which for smooth fields would only give exponential 
corrections, but for rough fields is part of the renormalization procedure) and 
the measure for integration over the fields needs to be converted from (com- 
pact) link variables to vector potentials. The main advantage of the lattice 
formulation is of course its intrinsic gauge invariance. In the continuum, even 
in a finite volume, the integral over the irrelevant gauge degrees of variables is 
hard to define unambiguously. We note that modifications in the integration 
measure can, in principle, be absorbed in the action. 

To compensate for the irrelevant higher dimensional operators in the action 
one can follow two routes. One is using Wilson's renormalization group, intro- 
ducing blocking transformations which for finite couplings and lattice spacings 
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define a xcnormalized trajectory along which the physical quantities remain 
constant Eil (giving a perfect lattice action). Its limit for zero coupling is the 
classically perfect lattice action E 2 ]. In a sense this can be seen, at least in prin- 
ciple, as a tree-level improved action to all orders in the lattice spacing, but 
through the blocking transformations is supposed to also correct for the other 
errors we mentioned abova. It may be expected, as demonstrated to one-loop 
order in the 0(3) model E3, that power-like cut-off effects in loop corrections 
are absent too. In this should lie the strength of classically perfect lattice 
actions, despite the need for truncations in actual numerical calculations. A 
large redundancy in the way one parametrizes these actions is reflected in the 
freedom of choosing the blocking transformations. One will try to keep the 
action as local as possible. 

The large redundancy is also clear in the Symanzik improvement scheme, 
which is based on perturbative computations of physical quantities. At tree- 
level there are many ways of choosing the action in terms of the gauge invariant 
Wilson loops that extend further than a single plaquette, such that the order 
a correction has a vanishing coefficient. The class of pactions we will consider 
is defined by (<> implies averaging over orientations Ej) 



S({c t }) = ^Tr{ Co (i-Q) + 2 Cl 



i i- 



3 C2 < 1 ~ 




4c, 



c 4 ( i- 



>} 



- y (co + 8ci + 8c 2 + 16c 3 + 16c 4 ) J2 ^{ F lu{x)) 
a 6 

+ — (c + 20 Cl - 4c 2 + 4c 3 + 64c 4 ) ^(P^F^x)) 2 

X.fl.V 

+« 6 T E Tr((^F„ A ) 2 )+0(a 8 ) • 



(3) 



We have added a 2 x 2 plaquette to the parametrization of the lattice action 
employed in ref. 8. (Note that sometimes in the literature conventions are used 
where ci and C3 are interchanged, e.g. refs. 10,14). We will motivate further on 
the advantage of adding this new plaquette. Using this action in perturbative 
calculations, one finds that the coefficients c» need to be corrected at higher 
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order in g^. This is achieved by computing a physical quantity to the required 
order and impose the vanishing of the 0(a 2 ) term to the relevant order in 
the coupling. This is called on-shell improvement and should be independent 
of the set of quantities required to remove the C(a 2 ) terms. Note that some 
redundancies appear due the invariance under field redefinitions. As was shown 
by Luscher and WeiszEj this allows one to choose C3 = 0. 

1.2 Tadpole improvement 

The most urgent question would be how small the coupling constant (and 
hence the lattice spacing) has to be chosen to be able to truncate the pertur- 
bative construction of these improved actions. Relevant for questions related 
to truncating a perturbative series is which definition of the coupling constant 
to take. Formally, lattice perturbation theory is defined in the bare coupling 
constant. Usually however, it helps greatly to convert the coupling constant 
to one that is defined at the scale of the process for which one is computing 
the perturbative corrections. This works well in the computation of physical 
quantities. However, one may doubt if this is valid for the computations of 
the improved action, as the corrections are by definition not physical. Still, a 
coupling at the scale of the cutoff would presumably perform better than the 
bare coupling. It has become customary to choose for this Parisi's mean field 
couplinglia, extracted from the plaquette expectation value, 

g 2 P = -aid, N) log(«o), u =< ^ReTrU(P) >^ 4 (4) 

where a is a constant, defined such that gp — go + 0(g$). It is clear that this 
is just one possible choice. Nevertheless in terms of this coupling, deviations 
from asymptotic scaling seem to be much smaller than in terms of the bare 
lattice coupling. One usually relates this coupling to one that is obtained Jay 
a resummation of tadpole diagrams, although this is hard to make precise 113. 
Similarly one may attempt to apply "mean field" type corrections to the im- 
provement coefficients that take into account the apparent large renormaliza- 
tion factors for the links, implied by the large difference between gp and go at 
moderate and large couplings. This has been known as tadpole improvement 
and is observed to work very well, in particular when considering the rota- 
tional invariance of the heavy quark potentialta. It should be emphasized that 
tadpole improvement can be seen as a rearrangement of perturbation theory. 
When only tree-level improvement coefficients are known the coefficients Ci are 
divided by a factor uq for each extra link (compared to the single plaquette). 

ci -> cx/uq, c 2 -> c 2 /uq, c 3 -> c 3 /wo, c 4 -> c 4 /uq ( 5 ) 
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For the only case sq_Ear where the one-loop improvement corrections to Cj are 
known, one adjusts cB the coefficients according to the expansion of gp in the 
bare coupling go. The implicit assumption of tadpole improvement is that 
corrections of O(gp) can be neglected. Although in many cases tadpole im- 
provement is shown to reduce the size of the perturbatiiffi corrections and to 
improve numerical results for very coarse lattices indeedliij, it remains a rather 
ad hoc procedure. As long as one needs to show the validity by comparing 
with the carefully extrapolated Wilson action results, one should shed doubt 
on the usefulness of this procedure in particular when one starts to "fudge" 
with alternate definitions of uq (for example extracted from more extended 
Wilson loops). What would be required is a well motivated definition of (a 
class of) improved lattice actions where similar to the Wilson action careful 
scaling studies are performed, that allow one to extrapolate results indepen- 
dently to the continuum and infinite volume. Working only at the coarsest 
lattices will probably prove to be insufficient to obtain reliable results. In 
particular systematic studies that use finite-size techniques as developed by 
the alpha collaboration-Ed for determining the running coupling constant and 
renormalization factorsEil seem mandatory for a careful study of the systematic 
errors involved. 

1.3 The setting 

Motivated by all these problems we set out to test improved actions in the 
context of finite volume spectroscopy. The clear advantage is that for volumes 
belowJ1.7 fm, results for the low-lying spectrum are know more or less analyt- 
ically Ej, both for the continuum and for the Wilson action. In such a case we 
can be precise in saying how much improvement is achieved. We can not reach 
the coarsest lattices employed so farEi Our volume has to be at least three 
lattice spacings across, which means that below 0.7 fm we can only reliably 
reach lattice spacings of the order of 0.2 fm. The data presented here will be 
for much smaller lattice spacing, where nevertheless scaling violations in mass 
ratios are particularly big. 

In order to also study the results of the improved actions analytically we 
introduced a new improved lattice action for which the gauge field propagator 
simplifies greatly. The explicit form of the propagator also will make clear that 
there are violations of unitarity (poles in the propagator with negative residue) 
at the scale of the cutoff. In the continuum limit these are harmless, in the 
same way as Pauli-Villars regulator fields are harmless below the scale_pf the 
cutoff. It does, however, cause difficulties in extracting a HamiltonianEl 
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2 The new improved action 

To study perturbation theory on the lattice we write U„{%) = exp(g„(a;)) and 
expand the lattice action to second order in q. At tree level we have, as for 
the Liischer-Weisz Symanzik action, C2 = C3 = 0, cf. Eq. 3. To remove the 
remaining 0(a 2 ) irrelevant operator we have the three constants Co, c\ and c 4 
at our disposal. One is eliminated by requiring the appropriate continuum limit 
without the need of rescaling the coupling constant. For the LW Symanzik 
action, corresponding to c 4 = 0, one finds cq = 5/3 and ci = —1/12. The 
disadvantage of this choice in perturbation theory becomes clear when writing 
the quadratic part of the action 

S 2 = ^-iTr[co(9^(x)-5^(a;)) 2 + 2 Cl {{2 + d^(d^(x)-d^(x))} 2 

+c 4 {(2 + d v ){2 + d^quix) - d.q^x))} 2 ] , (6) 

where the lattice derivative is the difference operator 

dp(p(x) = tp(x + fi) - ip(x) . (7) 

There is no "covariant" gauge condition that will make the gauge field prop- 
agator diagonal in the space-time indices. In analytic calculations this makes 
perturbative computations cumbersome, in_Darticular in the presence of a back- 
ground field, but of course not impossible^. Nevertheless, if we choose 

c 4 • c = c\ , (8) 

a condition invariant under the tadpole modification in Eq. 5, we can rearrange 
Eq. 6 by completing squares, from which we read off a simple gauge condition 

S 2 = - ^TcoTr^Or) (l + *(2+0+)(2 + M )) (l+z(2+dt)(2+d„)) d^q v {x)] 
+ E T </( a; )> = V^Y, d t (l+z{2+dl)(2+a it ))q M (xU9) 

where z = c\/cq. At tree level we find Co = 16/9, c\ = —1/9, c 2 = C3 = and 
c 4 = 1/144, whereas for ii / 1 we have z = — 1/16uq and Co = 1/(1 + 4z) 2 . We 
propose to call this new action the square Symanzik action, or square action 
for short. We find the following propagators 

Ghost : P(k) = ' 



ccTEa (4sin 2 (fc A /2)+4zsin 2 k x ) 
7^(1 + 4z cos 2 (fc M /2)) 



Vector : P^(fc) = ^ ^ . (10) 
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This simple form of the propagator allows us to demonstrate a general feature 
of improved actions, which have couplings in the time directions that are not 
nearest neighbor. We introduce the standard lattice propagators 



4sin 2 (fc /2)+ W 2 ' ± 2z ^ y (l + 4z) 2 y 

3 

loI = _i±i£ , w 2 = 4^sin 2 (fc J /2)(l + 4zcos 2 (fc i /2)), (11) 



which have a single pole in the fco Brillouin zone with the standard residue as 
would have been obtained from the Wilson action (except that u) = lo{z = Q) 
is replaced by u> s ). This allows us to write 



P = Z(P_ P + ), P„ = - ; )__J U ,^ P 



(l + 4g) 
(1 + 4zcos 2 (fc i /2)) i ' 



(l + 4z)(^P_-^P + + ^Po), (12) 



UJ 2 ' 



where the Z factor is given by Z = 1/ \J\ + 4zo> 2 / (1 + 4z) 2 , showing that 
there are unphysical propagating negative (and positive) norm states. This 
is a direct consequence of the fact that the transfer matrix, as defined in 
ref. 23, is not hermitian. The unphysical masses involved are at the scale of 
the cutoff (in lattice units 0(1)). For k = and uq = 1 one finds w_ = 
and u>+ = ujq — yl2- Although the unphysical states are interacting in a 
complicated way, and do not just appear in closed loops (thereby making it 
somewhat misleading to call them ghosts) they can be shown not to leadJxi 
imaginary eigenvalues and violations of unitarity in the low-lying spectrumcJ. 
In particular the low-lying spectrum for gauge theories in an intermediate 
volume is obtained by integrating out the non-zero momentum modes in a 
loop expansion and solving the resulting effective theory using a Raylcigh-Ritz 
analysis which includes imposing proper boundary conditions in field space 
implied by the gauge invariance E3. The effect of the spurious poles is integrated 
out at the same level of accuracy as the high momentum modes. Here we will 
only use our square Symanzik action to calculate the effective potential for an 
abelian constant background field, which greatly simplifies due to the diagonal 
structure of the gauge field propagator. Also it will be used to compute the 
Lambda parameter. 



7 



3 The effective potential 

The effective potential for a constant abclian background field C M is particu- 
larly simple to calculate since the background field will give rise to a shifted 
momentum, kn — > k„ + sC^/N^, where N„ is the size of the lattice in each 
of the four directions. For the modes q that commute with the background 
field one has s = 0, whereas for the two charged components with respect to 
the abelian direction defined by the background field one takes s = ±1. To 
one-loop order the effective potential is simply obtained by summing over the 
logarithm of the eigenvalues of the quadratic fluctuation operator, properly 
corrected for by the ghost contributions. We take Co = and N{ — N, such 
that 

v ^ c ) = w E lo gM fc + sC i N ) - lo g v(* + sC i N )} ( 13 ) 

fe,S (U 

The eigenvalues A(fc) are directly read off from the explicit expressions for the 
propagators in Eq. 10. 

Xgh(k) = v ^ "^4sin 2 (fc,/2)(l + 4zcos 2 (fc,/2)), 

K{k) = v ^ "(l + 4zcos 2 (fc M /2))A s , i (fc). (14) 

This formula holds equally well for the Wilson action (where z = 0). As fol- 
lows from the factorization of the propagator we can write X g h as the product 
of two eigenvalues. They will come in complex conjugate pairs when the mo- 
mentum gets too close to the edge of the Brillouin zone. Nevertheless splitting 
log Xgh in these separate contributions, that depend on fc as in the background 
field calculation for the Wilson action, allows one to perform the sum over fc 
exactly E3. In particular for Nq — > oo one finds the following explicit result 



V 1 {C) = iV^|4asmh ( 2u Jl+Az+Y +w^1+^ j +]Tlog(^)}, (15) 

up to an irrelevant overall constant. Here = l + 4z cos 2 [(2Trni + Ci) /(2N)] 
and oj = Lj(k = (irn+ ±C)/N) as defined in Eq. 11. For z — > this gives the 
well known result obtained for the Wilson action, whereas N — > oo recovers 
the result for the continuum theory lj. In figure 1 we compare the results 
of V\(C) for the square Symanzik action (lower two curves for N — 3 and 
N = 4), with those for the Wilson action (upper three curves for N — 3, 4 and 
6). The full curve gives the result of V\(C) for the continuum (we have chosen 
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Figure 1: The effective potential for a constant Abelian background field A\ = ^iCa^/N . 
The full line represents the continuum result. The lower two dashed curves are for the square 
Symanzik action with N = 3 and 4 (N = 6 is indistinguishable from the continuum curve). 
The upper three dotted curves are for the Wilson action with iV = 3, 4 and 6. 

C = (C, 0,0)). For the new improved action the result for N = 6 can already 
not be distinguished from that in the continuum at the scale of this figure. One 
might even fear that choosing uq ^ 1 will make the agreement worse. However, 
the effective potential is not a spectral quantity and deviations of V\(C) from 
the continuum can in principle be compensated by other corrections in the 
effective Hamiltonian for the zero-momentum modes. 

4 The Lambda parameter 

One can now follow the same strategy as in calculating the zero-momentum 
effective Hamiltonian for the Wilson action. The difficulty lies in convert- 
ing a non-hermitian transfer matrix, defined by the euclidean path integral, 
to a Hamiltonian. We postpone this to a future publication, where it will be 
shown that by a suitable field redefinition one can map the effective action with 
next-to-nearest neighbor couplings in the time direction to one with nearest 
neighbor coupling. After this transformation (whose Jacobian will give rise 
to corrections of odd powers in the lattice spacing), one can convert the path 
integral through a hermitian transfer matrix to a Hamiltoniano (this conver- 
sion will introduce corrections that are of even powers in the lattice spacing). 
At first, finding scaling violations that are of odd powers in the lattice spacing 
may seem rather puzzling, but it turns out they are required to exactly can- 
cel similar scaling violations that appear in some coefficients of the effective 
actionem 
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However, to determine the Lambda parameter for the square Symanzik 
action, it is sufficient to determine the effective action and study the finite dif- 
ference in the renormalization as compared to (say) the Wilson action in the 
limit of zero lattice spacing. Although in the renormalized action the tree level 
kinetic and potential terms, ^(d t cf) 2 and ^(F^) 2 , will have coefficients that 
differ by finite renormalizations (due to the breaking of Lorentz invariance), 
the difference of these finite renormalizations between the different regulariza- 
tions is expected to be the same for both terms. It should be noted that in 
principle one can perform a rescaling of the fields, but this would upset the 
periodicity (a consequence of gauge invariance) of the effective theory along 
the abelian constant potential. Indeed in background field calculations no field 
renormalizations appear. 

Consequently two numerically independent determinations of the Lambda 
parameter canJae extracted from this background field calculation, as for the 
Wilson action 123. In addition one can follow the alternative determination 
of the Lambda parameter through the heavy quark potential li-i All three 
determinations gave to a high degree of accuracy the same result for SU(2). 
The result for SU(3) was obtained with the latter method only. 



The index S 2 stands for the square Symanzik and W for the Wilson action. 
Details of both methods to compute these ratios will be presented elsewhere. 

5 Monte Carlo data 

In figure 2 we present the SU(2) Monte Carlo results for mass ratios in a small 
volume. Full details of the analysis will be presented elsewhere. We chose 
the volume such that the lattice artefacts in the mass ratios, the difference 
between the full curves (continmim) and the dotted curves (Wilson action for 
a 4 3 x oo lattice), were maximalEl For the Wilson action, at (3 = 4/g 2 = 3 this 
corresponds to a lattice spacing of approximately 0.018 fm. Another reason 
to pick these parameters was tp_.compare with earlier high precision data for 
the Wilson action by Michael E£l (crosses in the figure) , so as to make sure 
no errors were made in the implementation and in the measurements of the 
masses. We also improved here somewhat on the statistical errors. The data 
corresponding to the LW Symanzik action is represented by the triangles and 
for our new square Symanzik action by the squares. In both cases we used tree- 
level improvement only. The improvement is significant. For the LW Symanzik 
action the data is within two sigma of the continuum values. The results seem 



A S 2/A W 
A s2 /A w 



4.0919901(1) for SU(2), 
5.2089503(1) for SU(3). 



(16) 
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4 3 xl28 



□ square (3 = 3.3013 
a Symanzik (3 = 2. 374 
xo Wilson (9 = 3. 




Figure 2: SU(2) Monte Carlo data for mass ratios in a small volume on a lattice of size 
4 3 X 128, using the Wilson action (crosses from existing data of Michael), the LW Symanzik 
action and our new square Symanzik action. The lines give the analytic results: full for the 
continuum and dotted for the standard lattice action (N = 4). 



to indicate that the square Symanzik action is somewhat less effective, although 
the difference is not significant. A comparison at coarser lattices will be more 
interesting. The lattice results for the new action can be shown to confirm 
the value of the Lambda parameter computed from perturbation theory. The 
value of (3 for the Monte Carlo analysis with the square Symanzik action was 
chosen before the the perturbative calculations were completed. 

The mass ratios shown are all with respect to the scalar mass in the singlet 
representation of the cubic group. On the left we plot the square root of 
the finite volume "string tension" for one, two and three units of electric flux, 
whereas on the right we consider the tensor doublet E + and triplet , which 
become degenerate in large volumes when rotational invariance is restored. 
The finite volume "string tension" is simply defined by dividing the energy of 
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the electric flux states (also called torelon masses by the length of the box.. 
The electric flux is defined with the help of twisted boundary conditions E3. 
The Wilson loop operator that closes due to the periodic boundary conditions 
can be seen to create electric flux. Two (or three) units of electric flux are 
obtained from loops that wind around two (or three) directions of the torus. 

6 Outlook 

Our contribution has been modest as we have not probed very coarse lattices. 
The Monte Carlo data presented here were for a lattice spacing of approxi- 
mately 0.02 fm. New Monte Carlo data will be presented elsewhere, studying 
lattice spacings of approximately 0.12 fm. We would also like to stress that we 
used mass ratios to test improvement. It has been well know that the Parisi 
mean field coupling constant improves the approach to asymptotic scaling quite 
well. We see this as a separate issue, quite (but not completely) unrelated to 
the issue of scaling. The reason is that the scale in lattice calculations is usu- 
ally set by one of the masses (in pure gauge theories by the string tension) 
anyhow. 

From the practical point of view we presented an alternative action, moti- 
vated by the requirement to simplify the perturbative calculations, rather than 
to minimize the lattice artefacts. Nevertheless, comparisons will give a clue on 
how big the higher order corrections are. Of course at tree level it is trivial 
to write down many types of improved actions, so we have set out to compute 
the one-loop improvement coefficients too, to bring our new square Symanzik 
improved action to the same footing as the Liischer-Weisz Symanzik improved 
action. These results will be presented elsewhere. This allows one to test the 
"universality" of tadpole improvement. Again we wish to stress that the most 
interesting tests are those for the scaling of mass ratios. 
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